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ABSTRACT 

The article includes the results of theoretical studies of the accuracy of geodetic height survey and 
marks points on the Earth's surface using satellite technology. The dependence of the average square 
error of geodetic heights difference survey from the distance to the base point was detected. It is being 
proved that by using satellite technology one can define vertical displacements of the Earth's surface 
with accuracy sufficient for practical application. It was revealed that the use of positioning methods 
after the signals of global navigation satellite systems GPS and GLONASS (GNSS) for geodetic purposes 
has begun in Russia in the beginning of 90-s of the last century. Their essential advantages compared to 
traditional survey methods were revealed. These include a wide range of accuracies (from tens of meters 
to millimeters), independence from weather, time of day and year, absence of necessity for mutual 
visibility between points, high automation and, consequently, efficiency, ability to work continuously 
and in motion. These qualities have led to high performance and economic efficiency of GNSS, which is 
particularly noticeable in remote and unpopulated areas occupying a large part of our country. 
Currently, according to the results of the measurements by double-frequency satellite receivers the 
mean square error of coordinate increment calculation was reached equaling to 3 mm+1-10''bD, where 
D is the distance between the satellite receivers. But exceedances between the very same points can 
be obtained with the mean square error of 10-30 mm, which greatly increases along with the increase 
of the distance (D) to tens of kilometers. In traditional geodesy height and exceedances surveys are 
conducted relative to the surface of quasigeoid, which means that measurements are based on physical 
principle of measurements. As a result, geodetic networks built by traditional methods can be divided 
into planimetric (B and L) and vertical control HI networks, which are almost unrelated. In this regard, 
for the purpose of enhancing the effectiveness of geodetic application of the satellite leveling, 
methodology and technological support of satellite measurements require improvement. 
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Introduction 

Numerous effects on the Earth's surface occurring during aboveground and 
underground construction and production of liquid and solid commercial minerals, 
lead to its deformation (Koneshov et al., 2016; Bol’shakov, 1997; Drobyshev, 2006; 
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Nepoklonov, 2010; Schultz & Winokur, 1969; Nesenyuk et al., 1980; Dmitriev, 1997; 
Timochkin, 2013). Timely and expeditious detection of reliable deformation 
parameters of the Earth's surface allows taking effective measures for engineering 
protection of the objects located thereon (Klyushin & Kravchuk, 2009; Gairabekov et 
al., 2009; Gairabekov, 2011a; 2011b; Gairabekov & Pimshin, 2010; 2011). Information 
on deformation processes is necessary for the proper design of new construction sites 
as well (Li & Jekeli, 2008; Peshekhonov et al., 1995; Fagerlund, 1977; Hirt et al., 2010; 
Kudrys, 2009; Tsodokova et al., 2014; Vasil’ev et al., 1991; Troitskii, 1994). For quick 
data updates of the relevant accuracy advanced geodetic measurement means should 
be used, which include geodetic satellite technologies (Semenov, 2012). 

Materials and Methods 

Satellite technologies allow accurate fixation of X and Y planimetric coordinates. 
Thus accuracy of coordinate increments with the use of double-frequency satellite 
receivers equals to 3 mm+l-10-6D, where D is the distance between station points of 
satellite receivers. Vertical control position of points is considered to be fixed with 
much less accuracy, with the mean square error of 10-30 mm. However, this problem 
is still under-explored with both theoretical and experimental study of this problem 
remaining the relevant objective. We are going to undertake theoretical study of the 
possibility of using satellite technologies to detect vertical displacements of the Earth's 
surface. 

The article (Klyushin & Kravchuk, 2009) gives the derived formula, which allows 
computing the geodetic height of a point based on the results of satellite 
measurements: 

H = [X 2 + Y 2 + Z 2 ( 1 - a 1 R 2 r~ 3 )~ 2 (l + b^r' 3 ) 2 ] 1 ^ 

e 2 Z 2 R~ 2 

~ a f “ Z 2 R~ 2 + (1 - a 1 R 2 r~ 3 ) 2 (l + b 1 Z 2 r~ 3 )~ 2 

where % = ae 2 (l — e 2 ) 3 /z and b 1 = be' 2 ; a and b are ellipsoid semiaxises adopted for 
the processing of geodetic measurements and establishing a geodetic coordinate 
system; e 2 and e' 2 are values of eccentricities; 

R = VXMW 2 

r = 7(1 -e 2 )(Y 2 + F 2 )+Z 2 


( 1 ) 


Formula (1) can be represented as follows: 

H = P — N, (2) 

where N is the curvature of the first vertical; 

P = 7V 2 +Y 2 +Z 2 ( 1 - a 1 I? 2 r- 3 )- 2 (l + b^r- 3 ) 2 


If we take the differential equation (2) we get: 

dP dN dP dN dP 
dH = —dX -—dX + —dY -—dY - —dZ 
dX dX dY dY dZ 

The following designations are accepted: 

L = 1 — a 1 P 2 r -3 , (5) 
t 2 = 1 + b 1 Z 2 r~ 3 , (6) 

t 3 =^f,(7), (7) 

t 4 = S + p (8), (8) 


dN 

dz‘ 


-dZ 


(3) 

(4) 
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Vertical curvature we present in the following way: 

a 

N = 

£3 

U 


Formula (4) can be represented as follows: 


P= X 2 + Y 2 + Z 2 t~ 2 t 2 


cLH = 


After some transformation of the formula (4) we get: 

2 ^ V7* or, + 2 v 72 3fli (i - e 2 )tlXZ 2 R 2 


X 3b 1 (l - e 2 )t 2 XZ 4 2a 1 t^XZ 2 
P Pt^r 5 + Ptfr 3 

at 3 XZ 2 


Pt^r 5 


+ 


ae 2 XZ 2 


"2iCLCL-\ £1 toX 


t 4 / 3 (t 4 - t 3 ) 3 /2R4 tf t 4 /3 (t 4 - t 3 ) 3 /2r 3 


+ ■ 


.1/3 
L 4 

3aa 1 t 1 t 3 (l — e 2 )XR 


3 / 

2 A VD2 


t;-(t 4 -t 3 ) /2R* 


L 4 L 3 

’ll tits 


tft 4 / 3 (t 4 


t 3 ) 3 /: 


27 ° 


2 

dA 


tft 4 / 3 (t 4 - t 3 ) 3 // 2 r b 


+ 


F 3^(1 - e 2 )t 2 FZ 4 Za^iFZ 2 3^(1 - e 2 )t 2 2 FZ 2 F 2 
P Pt?r 5 + Pt, 3 r 3 Pt?r 5 


+ ■ 


ae 2 YZ 2 


2aa 1 t 1 t 3 Y 


M2, 


t 3 ) V 2 P 4 


+ 


+ 


(t 

c 4 \ l 4 u 3 J ~ 11 ^4 

361 % *^£ 3(1 — e 2 )YR 2 


ty 3 a 4 


t 3 ) 3/ 2P 4 
3 ab 4 tl t 3 YZ 2 


t|t 4 ' (t 4 


t, t 4 / 3 (t 4 - t 3 ) 3 /: 


u 2 ^4 

Zt| 
Pt 2 ' 


l 4 

3b 1 t 2 Z 

Ptfr 5 


3 / /2 r b 
5 


t|t 4 /3 (t 4 - t 3 ) 3 /27 
ae 2 Z 


1-2 *- 4/3 ( 

dF 


t 3 )'■ 


Ptfr 3 


.1/3 


(t 4 


3am t-, t 3 ZJ ? 2 


t|t 4 / 3 (t 4 — t 3 ) 3 / 2 r ; 


■ + 


3 ab 4 tf t 3 Z 3 


t|t 4 / 3 (t 4 — t 3 ) 3 // 2 r 


t 3 ) 3/ 2 F 2 

dZ 


+ ■ 


2^3 


atoZ 


t 4 4/3 (t 4 

(9) 


t 3 ) * 2 F 2 


Since there is a small quantity of difference between t^and t 2 coefficients and 


the unit, formula member (9) can be represented as follows: 


Z(l-a 1 fi 2 r -3 ) 


2a,ZR z 2b 1 Z s 


( 10 ) 


P 2 !^ P(l+b 1 Z 2 r -3 ) 2 P Pr 3 Pr 3 

And given the (10) we present the formula (9) as follows: 

dH = (£ + A)d* + (£ + fl)dF + (^+c)dZ, (11) 


where 


4 


3^(1 - e 2 )t 2 ZZ 4 2a x t\XZ 2 
Pt 2 r 5 + Ft 3 r 3 


at 3 XZ 2 

t 4 4 / 3 (t 4 - t 3 ) 3 / 2 F 4 
3ab-,tl t 3 XZ 2 
t|t 4 / 3 (t 4 — t 3 ) 3 // 2 r 5 


3a 4 (l - e 2 ')t\XZ 2 R 2 ae 2 XZ 2 

Pt 3 T 3 + t l/3 ( - t4 _ t3) 3/ 2/?4 

2aa 1 t 1 t 3 X 3aa 1 t 1 t 3 (l — e 2 )XR 2 

t|t 4 / 3 (t 4 — t 3 ) 3 ^ 2 r 3 tft 4 / 3 (t 4 — t 3 ) 3 / 2 r 5 

; (12) 
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B = — 


3^(1 - e 2 )t 2 yz 4 2a 1 t%YZ 2 


Ptfr 3 


Pt^r 3 


3a x (l - e 2 )tlYZ 2 R 2 
Pt\r 3 


+ 


ae 2 YZ 2 


aUYZ 2 


2 am t-, t 3 Z 


t 4 /3 (t 4 - t 3 ) 3 /2/? 4 tft 4 /3 (t 4 - t 3 ) '2r 3 

; (13) 


+ ■ 


1/3 , 

4 

3aa 1 t 1 t 3 (l — e 2 )YR 2 


- t 3 ) 3/2 i? 4 


L 4 L 3z 

V 2 t 3J 


t| t^ 3 (t 4 


ts) 3 /: 


27° 


t|tl / 3 (t 4 - t 3 ) 3 // 2r b 


C = - 


2a t ZR 2 

Pr 3 


3b,UZ s 3a,R 2 Z 3 


Pr 3 


Pt 2 r 5 


Pt?r 3 


3aa-, U t 3 ZR 2 


t| t 4 /3 ( f 4 ~ t 3 ) 3 / 2 r ; 


■ + 


ae 2 Z 


t* /3 (t 4 - t 3 ' 3 
3ab 4 tl t 3 Z 3 


t 3 ) 3/ 2R-. 


■ + ■ 


.4/3 


citoZ 


(t 4 - t 3 ) 3 /2i?2 


t | t 4 /3 ( t 4 ~ ts) 3 ^ 27 "" 


(14) 


Next we can move from differentials to finite increments and then to mean square 
errors, obtaining the following from (11): 

= + m3 + (^ + 5 ) m r + (j;+c) m|(15) 


The analysis of accuracy of point coordinates survey using modern satellite 


technologies reveals that generally TTLx ~ 777y — 772^ — 772^. Therefore, the 
formula of accuracy assessment of geodetic height calculation of a point can be 
presented as follows: 


m H = m K 


y 


X 2 + Y 2 + Z 2 24Z 25Z 2CZ 

-^+ -r- + i4 2 + B 2 + C 2 (16) 

P 2 P P P y J 


Results and Discussion 

The performed formula analysis (16) showed that with an error not exceeding 4%, 
it can be written as: 

m H = m K 


X^+Y^+Z^ 

Since, ——-~ 1 the remaining members are small. 

By way of example, we get the following values in the middle line: 
B = 45°; L - 37; H = 200m 

X = 3608020 m; Y = 2718839 m; Z = 4487489 m 


X 2 + Y 2 + Z 2 


P 2 


= 1.023 


2 AX 
~P~ 

2 BY 
— 

2CZ 

— 


= -0,0048 
= 0,0089 
= 0,0397 


i4 2 = 0,000018; K = 1,0237 
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B 2 = 0,000108 
C 2 = 0,000791 

At the equator we get; B — 0 °;L — 37°; H — 200m 

X = 5093965 m; Y = 3838578 m; Z = 0 m; 


X 2 + Y 2 +Z z 


P 2 


2AX 

— 

2 BY 
~P~ 

2CZ 


= 0 


= 0 


- o 


■= 1 


i4 2 = 0 
B 2 = 0 
C 2 = 0 


In the area of the Arctic covered settings will have values: B — 72°; L — 37°; H = 
200 m 

X = 1578909 m; Y = 1189793 m; Z = 6043875 m; 


X 2 + Y 2 +Z 2 
P 2 


■= 1,0012 


2 BY 

— = 0,0036; K = 1,0593 
2CZ 

— = 0,1163 

i4 2 = 0,0060 
B 2 = 0,00009 
C 2 = 0,0037 

It should be noted that the geodetic heights are tied to point coordinates identified 
by means of satellite technologies and the chart of complex and non-linear form. 
However under the equally accurate Cartesian coordinates, mean square error of 
geodetic height calculations is essentially independent of a point position. 

To assess the accuracy of the results of satellite measurements, considering 
formula (2), we state: 

h = h 2 -h 1 =p 2 -p 1 + n 1 -n 2 


(17) 
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where 


P^ = 


, , . (1 + \ Zfr- 3 ) 2 

2 i V2 i 72 y_ 111 _' 


X'l + Y 2 + Z 1 


(1 — a-tRlr-^ 3 ) 2 


P 2 = 


(1 + b x Z\r 2 3 ) 


v2 I yT. I 72 v ~ ' 2 2 

2 2 2 (1 - aii?|r 2 - 3 ) 2 


N, =■ 


1 - 


e 2 Z 2 fl~ 2 


r2 n-2 1 (1 a-iRfvi 3 ) 2 


zfflr 2 + 


(1 + MiV 3 ) 2 


iV 9 = 


e 2 Z 2 /? 2 2 


Z 2 f?7 2 + 


(1 — a t Rlr 2 3 ) 2 
(1 + ^Z|r 2 - 3 ) 2 


Here below we give the expression of a total differential of geodetic heights 
difference: 

dh = (| + 21^ dX, + (| + 1?!) d^ + (| + C,) dZ ! - (| + i4 2 ) dZ 2 - (| + B 2 ) dY 2 - 

(|+C 2 )dZ 2 (18) 

The coordinates of the second point (Z 2 , b 2 , Z 2 ) can be defined using the 
coordinates of the first point (Z 1( Y 1 ,Z 1 ) and coordinate increments (AX, AY, AZ) phase 
measurements. Accordingly, 

dX 2 = dX 1 + dAX) dY 2 = dY 1 + dAY; dZ 2 = dZ 1 + dAZ, (19) 

dH = (^-^ + A 1 -A 2 )dX 1 + (^- Y f 2 +B 1 -B 2 )dY 1 + (^- Z f 2 +C 1 -C 2 )dZ 1 - 
g + A 2 ) dAX - g + fl 2 ) dAb - g + C 2 ) dAZ. (20) 

Next we go from differentials to finite increments, and to mean square errors 
afterwards, so that we get the following from (20): 


= (S“S +i4i_i4z ) m ^ + (!~J +Si-S2 ) m ^ + (t~S +Ci - c 2 ) 

g + ^z) m Ax + g + #2) m Ar + g ~ ^ 2 ) m AZ (21) 


Subject to equally accurate positioning of point 1 m Xi = — m Zi = ml 

and coordinate increments m\ x — m\ Y = m\ z = m\ formula (21) is as follows: 
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m h 


[(l-l) 2 - S-S 2 - (l-S 2 - 2 (Ai-A 2 ) (l-l) + 2( Bl -B 2 ) (|-g 4 
2(C 1 -C 2 ) (|-|) 4- (A^A,) 2 4- (B 1 -B 2 ) 2 (C 1 -C 2 ) 2 ] m 2 + [g + A,)' + g + B,)' 

g + cglm 2 (22) 


4- 


Our formula (16) analysis shows that with high level of accuracy the value is 


GJ + Az) +(| + B 2 ) +(|+ C 2 ) 


1. 


(23) 


With this in mind, the formula (22) can lead to: 


mt 


(M)h (M)' + (M)' + 2 W , - , 2 > g -1)+2 (Bl - g-l) 

2(C, - C 2 ) (| - |) + (Ai - A 2 ) 2 + (B, - B 2 ) 2 (C 1 - C 2 ) 2 ] m 2 + : 


--- 4- 
ypi p 2 J 

(24) 


The analysis of the members constituting the formula (13) was performed under 
the condition that P 2 = Pi 4- A Pj which provided a minimum correction value of AP. 
After the appropriate conversions one can get the value AP equaling to: 

(25) 


_ X^X Y^Y ZiAZ z|(5 2 -5 3 ) AX 2 +AK 2 +AZ 2 


Pi 


Pi 


Pi 


Pi 


2Pi 


then the formula (24) is written as follows: 


m 


2 _ 


r ml + m\ + A„ 


(26) 


where 


A,„= 2(A, - A 2 ) g - |) + 2^ - B 2 ) (| - I) + 2(C, - C 2 ) g - |) + (A, - A 2 ) 2 + 
(fl, - B,) 2 ^ (C, - C 2 y + «* + + S 2 AP 2 , (27) 


Pi 


Pi 


Pi 


S 2 = AX 2 + AY 2 +AZ 2 


(28) 


The A,„ amount is negligibly small so we neglect it. For example, if B — 45°; L = 37; 


H = 200m A m = 0,000674 


Consequently, the expression of accuracy assessment of geodetic heights 
difference can be written as: 


m h 



(29) 


where ~ R 3 , R 3 is the mean radius of the Earth. 

On the basis of the accuracy analysis we conclude that at relatively large 
distances between defined points, base station coordinate error has a significant 
impact on the calculated values of the difference between the geodetic heights. The 
coordinates of the base station are calculated with the mean square error equaling 
TTl ^ =5 m, 5 m, and coordinates increments are calculated with the mean square error 
5 mm. The dependence of the accuracy of difference between geodetic heights 
calculations and differences of the coordinates is presented in Figure 1. The diagram 
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shows that with the distances from a base station not exceeding 5-6 km coordinate 
error has a significantly less effect on the accuracy of heights difference calculations 
than mean square error of coordinate increments calculation. 



Figure 1. Diagram of dependence of the mean square error of identifying the difference 
between geodetic heights distance from the distance to base station 

As it follows from the formula (18), the mean square error of identifying vertical 
displacements of points equals to the mean square error of identifying coordinate 
increments of this point. In the case we are considering it equals 5 mm. With the use 
of modern satellite double-frequency receivers of coordinate increment one can identify 
with the mean square error of 2-3 mm (Gairabekov et al., 2015a; 2015b), hence vertical 
displacements (precipitations) are calculated with the same accuracy. Such accuracy 
is sufficient to solve the above-mentioned problems: 

- calculations of vertical displacements of the Earth's surface of anthropogenic 
nature as a result of overground and underground construction; 

- calculations of vertical displacements of the Earth's surface during the 
extraction of liquid and solid mineral resources; 

- calculations of vertical offsets of seismotectonic nature. 

Conclusion 

Thus, the accuracy of calculation of vertical displacements by means of satellite 
technologies with distances between points up to 6 km is commensurate with the 
accuracy of calculation of horizontal displacements. Therefore, deformations 
(including vertical displacements) can calculated with use of satellite technologies with 
precision sufficient for practice in a relatively short period of time, which allows taking 
that prompt and effective measures for protection of engineering facilities located at 
the territories exposed to deformation processes, as well as obtaining information 
necessary for the design of new facilities. 
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